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This grant is to provide support for gathering and analysis of seismic, grav- 
ity and electrical resistivity data that are believed to be most relevant to 
development of earth'^uake premonitory models of the crust. These methods 
have been chosen to best complement the crustal distortion information meas- 
ured by the NASA extraterrestrial geodetic techniques, in particular, the ARIES 
project, Simultaneous measurement of these geophysical parameters with the 
NASA extraterrestrial geodetic technique, is considered essential to develop- 
ment of a physical model of the crustal processes preceding earthquakes. 

Seismic InvesLigations 

Preparing data bases in an easily accessible formal to researcljers is 
important in an area with a high seismic hazard such as Southern California. 
Hand-witten records of seismic data from the beginning of the Sei.smological 
Laboratory comprise the sole data base of seisinicity in Southern California 
since 1928, However, systematic re»="?arch on this data has been \drtually 
impossible because of the large task of putting the readings into a computer- 
compatible format. This ta.sk, involving the punching of seismic data, ha.s now 
been done for the most recent, and most voluminous, part of the data base. 

A card image computer tape is now available with earthquake epicenter 
and P- and S-phase data for local earthquakes in Southern California during 
the 10-ycar period of 1966 through 1975. The tape is readable on an IBM 
370/158 computer using standard FORTRAN READ statements. The tape is 
currently being used for several research projects including event relouallon, 
velocity, P-delay. and magnitude studies. This research re.source is being made 
available to all scientific investigators. 

The Seismological Laboratory was founded by the Carnegie Institute of 
Washington in 1927 and management was transferred to the California Institute 
of Technology (Caltech) in 1936. Although some seismographic stations of the 
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Caltech network were established as early as 1926, routine and systematic epi- 
central detenninations started only in 1932, with 7 stations reporting. By 1966, 
the beginning of the data tape, the network had expanded to 19 stations, By 
.1972, 39 stations were being read for earthquake phases. During the period 
from 1972 through 1975, the last year of the tape, a. cooperative program was 
established with the U. S. Geological Survey and a major expansion of the 
Southern California network increased the total number of stations to 103. 

Currently, the number of stations in Southern Calif ornie exceeds 150. This 
exponentially expanding data base has threatened to saturate even routine 
analysis of local earthquakes and led to the development of more automated 
methods of recording and analysis. These new methods include the Caltech 
Earthquake Detection and Recording (CEDAR) system, which has been In opera- 
tion since January 1, 1977. Automated earthquake analysis systems need u 
means to extrapolate their results to past data which, until the generation of 
this magnetic tape, have been inaccessible for computer processing, 

The magnetic tape contains approximately 150,000 card images represent- 
ing more than 200,000 individual readings, more than one-half are from the last 
2 years (1974 to 1975) of the 10-year interval, representing an increase in both 
the number of stations and the sensitivity of detection/location thresholds. 
The data tape is formatted chronologicall}' by event such that '■-ach ovent 
begins with a hypocenter description card followed by P- and S-phase card:; for 
that event. Some phase data follows blank hypocenter cards indicating tliat the 
phase data was deemed insufficient to attempt a hypocentral location. One 
yeeu- of data is linuted to one tape file, and there are as many records within 
each file as necessary: records are limited to 1000 cards or less. During the 
1966 to 1975 10-year span, hypocenters have been computed with difrerent 
location programs and different assumptions as to velocity models and station 
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delays. Hopefully, these variables have been modified ■with time tovrard 
improved hypocentral locations. But those studies depending on the uniformity 
and relative precision of hypocenter calculations during the 10-year span will 
probably require recomputation with uniform location parameters. Indeed, one 
of the major incentives for making the phase data available to computation is to 
remove the dependence of such studies on the particular method and assump- 
tions of h}'pocentral computation which changes over a long time span 
(Wliitcomb, 1978). 

The closest major fault to the NASA long-baseline geodetic sii.e Mount Otay, 
which is used by the laser-ranging SAFE and VLB! ARIES projects, is the Elsinore 
fault. The seismicity of this fault bears directly on the formation of tectonic 
models for distortion in this region, and not surprisingly on the seismic hazard 
in the San Diego area. A relatively small Ml = 4.8 earthquake and its aftershock 
series on the southern portion of the Elsinore Fault Zone in eastern San Diego 
County, California, provided a rare opportunity to study this area which has 
been subjected to variable tectonic interpretations in the past. Within IZ to 26 
hours after the main shock, a network of four portable seismograph stations 
was established around the main event near Aqua Caliente Springs to supple- 
ment the stations of the Southern California Seismographic Netwoik. Four days 
after the main shock, seven additional portable seismograph stations v.ern 
in.stalled. In addition to the main event, 45 subsequent events were studied, 
ranging in magnitude from about 1.0 to 3.7. Of these. 36 could be termed aft- 
ershocks by their close proximity to the main event, whose proper location was 
determined by analysis of the aftershock series, Of the two branches of the 
Elsinore Fault in this region, the south branch is associated with the earth- 
quake series. Focal mechani.sms are consistent with right-lateral strike-slip 
along the south branch, with northeast dip at latitude 32“51'N. These conclu- 
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sions are supported by hypooentral locations. Thrust activity on thr tvo fault 
branches may be developing a horsf between them, accounting for elevation 
and tilt change^,’ observed near Aqua Caliente (Allison et al., 1978), 

A great deal of potential information about an earthquake is contained 
■within the coda or the energy after the body phases recorded at a seismic sta- 
tion. Many questions relating properties of the coda to the size of the earth- 
quake and to the properties of the medium of the crust in the vicinity of the 
earthquake have been unanswered until recently. Research on these questions 
has produced a relationship between the seismic moment, Mo> of shallow local 
earthquakes, coda amplitudes, and the total duration of the signal, t, in 
seconds, measured from the earthquake origin time. The theory assumes that 
the end of the coda is composed of backscattering surface waves due to lateral 
heterogeneity in the shallow crust. Using the linear relationship bv^U, ten the 
logarithm of Mq and the local Richter magnitude Mi, obtained a relationship 
between Mi and t, of the form: Mi = o.q + ay log t + + / (0 , where Qq* 

ag are constants depending on an attenuation parameter (effective Q) and 
geometric spreading: and f{t) is a function of the scattering process. This 
relationship is different from the empirical one generally used 
Ml = tto + ay log r + aziloy r)^ + agA , where r is the duration measured from 
the first P arrival time and A is epicentral distance in kilometers. In the 
theoretical relationship, the dependence on epicentral distance is implicit in f. . 
The theoretical relationship is used to calculate a coda magnitude Me that is 
compared to Mi for Southern California earthquakes which occurred during the 
period from 1972 to 1975. This comparison is made independently at six sta- 
tions of the CIT network. At all stations, a good linear fit {Mi = Cq CyMc) is 
obtained. The standard errors range from 0.2 to 0.3 and the correlation conUi- 
cients from 0.80 to 0.90. Once station gain is accounted for, station correlation 


“ 6 ” 


ORIGINAL PAGE iS 
OF POOR QUALITY 


terms are less than 0.17 magnitude unit when comparing Mi and Mc> Me calcu- 
lation is not limited to a duration measurement but can utilize the entire earth- 
quake coda in order to increase by many times the statistical confidence in an 
estimate of an earthquake’s magnitude. Temporal variations in the parameter 
Q may be fundamentally related to models of the process of preparaticin for 
earthquakes (Suteau and Whitcomb, 1979). 

All researchers in tectonophysics that are trying to measure geonhysicai 
paramei.ers of crustal distortion have a need to be informed in a timely manner 
of tlje level of seismicity in their particular area of interest. In Southern Cali- 
fornia, the sole source of this information is the Seismological Laboratory, 
While the seismicity reports out of this institution have been more rapid than 
any other similar agency, there is usually a delay measured in years from the 
time of the earthquake to the time of the report detailing its occurrence. In 
the interest of impro^nng the timing of the seismicity reports, preliminary 
monthly seismicity reports have been instituted beginning in 1978 that have a 
delay that is now measured in months instead of years. This service continues 
at the present time. Ultimately, the goal is to publish reports for a month’s 
seismicity no later than one week from the end of the report period (Whitcomb 
et al., 197B). 

Gravity Investigations 

The Southern California gravity monitoring project, begun in May 1974, is 
intended to coordinate gravity measurements with the long- baseline three- 
dimensional geodetic measurements of the AIRES (Astronomical Radio Inler-- 
ferometric Earth Surveying) project which uses radio interferometry with 
extra-galactic radio sources. Gravity data from 28 stations, shown on Figure 1, 
are monitored on an approximately one- to two-month basis. Appendix I shows 
the gravity data from the first five years of observation of the networks. 
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Different symbols correspond to different meters used in the survey. The dou- 
ble line represents plus and minus two standard errors of the mean of a five- 
point weighted moving average. Several lines of evidence indicate that the four 
standard error interval is a good estimate of the actual gravity value which can 
be used to test the significance of a change in gravity. 

Electrical Resistivity Investigations by Means of MagnetoteUurics (MT) 

Magnetotelluric (MT) determinations of crusted electrical resistivity is one 
of the most promising techniques for determining the depth dependence of cru- 
stal distortion. It, therefore, has a direct bearing on the formulation cf distor- 
tion models. Early in the grant period, substantial modifications were made to 
the MT hardware which was developed under earlier grants. These included 
redesign of the electric field pre- amplifiers and major changes to the data- 
logger logic to enable the 9-track tapes to be processed more efficiently. 
Software was implemented on the departments NOVA and ECLIPSE computers 
which enabled high-speed plotting of both raw and digitally filtered data. Only 
two sites were occupied during the first six months of 1978. This was partially 
due to problems in the electric field preamplifiers which were not fully resolved 
until July, In addition, the China Flats and Lytle Creek sites were inaccessible 
for the fir.st six months of 197B due to heavy rainfall washing out the access 
roads. Good data was obtainved from the three MT sites during the months of 
July, August, September and November of 1978. Processing has been performed 
on this data to obtain maximum and minimum apparent resistivity curves ip the 
frequency range 0.001 - 30 Hz, Significant changes in high frequency apparent 
resistivity estimates have been observed at all three sites. These have been 
interpreted in terms of reduced near surface resistivities due to heavy rainfall, 

A study of the effects of magnetic source field polarization was performed 
on data collected from the three sites, It was discovered that highly pc’arized 
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source fields can cause a slight upward bias in the diagonal elements of the 
impedance tensor. Software has been developed to remove the biasing effect of 
noise on the resistivity estimates. The technique utilizes a singular valued 
. decomposition of the matrix of Fourier coefficients of the electric and, mag- 
netic fields. This technique considerably improves the repeatability of the 
resulting resistivity estimates. 


We have developed a technique which can be used to estimate magnetotcl- 
luric impedanct- from selectively polari?:ed source fields. The technique 
involves digital filtering of time-series dnla, followed by estimates of the ampli- 
tudes and phases of single cycles, utilizing a simple "maximum/minimum" 
approach. This yields impedance estimates that are at least as consistent as 
estimates obtained using the more usual Fast Fourier Transform/bandwndth 
averaging technique. When the source fields arc subdivided into intervals en 
the basis of polarization direction, the diagonal elements of the impedance ten- 
sor determined in each interval show an upward bias. We interpret this as due 
to instability in the solutions for the impedances which occurs when the source 
fields are highly polarized. In addition, one of the tensor elements shows a con- 
sistent upw'ard trend w'hen the source field is polarized close to the easi-west 
direction. Since we observe this at tliree separate sites, we suggest that it may 
be due to either a noise source which is polarized in the east-west direction, or 
to a lack of uniformity in the horizontal source fields which are polarized in this 
same direction (Lienert, 1979). 


Radon Investigations 

It is difficult to find other long-term geophysical measurements that arc 
believed to be relevant to the proce.ss of preparation of earthquakes with long- 
term stability comparable to that of NASA’s long-baseline geodetic techniques 
such as AIRES. One such technique is the monitoring of radon emanating from 
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the crust. This grant has been used to .support the initial effort to develop an 
automated radon monitoring system for installation near ARIES sites. The first 
such site is near the ARIES JPL site in Pasadena. Up to this time data have been 
gathered from 20 months of near real-time (3 samples per day) radon monitor- 
ing at this hard-rock site in the Transverse Ranges of Southern California. An 
annual cycle is evident in the data which is attributed to thermo-elastic sliains 
in the vicinity of the borehole; site. Between 1 April 1977 and 31 October 1973 
there vmre 11 earthquakes with magnitude 2.0 within 25 km of the monitor- 
ing site. Three of these events appeared to be preceded by precur'?ory .■signals, 
four were preceded by "possible" precursory signals, and four were not pre- 
ceded by any apparent precursors. Before the 4.6 M Malibu earthquake of 1 
January 1979, a possible precursory signal sequence of 40-50 days duration was 
observed (Shapiro et al., 1978). 

Other Work 

During the grant period, contributions were made to educational and 
scientific publications including 1) a popularized article on earthquake pr^dic'- 
lion (Whitcomb, 1978), 2) an invited paper in a National Academy of Scienres 
study on the impact of technology on geophysics (Whitcomb, 1979), 3) a contri- 
bution to NASA Technical Paper 1464, and 4) a study of real time data 
management/analysis systeni for operational earthquake prediction (Lehman 
et al„ 1979). 
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Appendix I. 

Gravity data from the first five years of operation of 
the Southern California Gravity Network. 
Station locations are shown in Figure 1, 
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